Fingering instabilities of a reactive micellar interface 
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We present an experimental study of the fingering patterns in a Hele-Shaw cell, occurring 
when a gel-like material forms at the interface between aqueous solutions of a cationic surfactant 
(cetyltrimethylammonium bromide) and an organic salt (salicylic acid) , two solutions known to form 
a highly elastic wormlike micellar fluid when mixed homogeneously. A variety of fingering instabil- 
ities are observed, depending on the velocity of the front (the injection rate), and on which fluid is 
injected into which. We have found a regime of non-confined stationary or wavy fingers for which 
width selection seems to occur without the presence of bounding walls, unlike the Saffman- Taylor 
experiment. Qualitatively, some of our observations share common mechanisms with instabilities of 
cooling lava flows or growing biofilms. 

PACS numbers: 83.60.Wc, 47.20.Gv, 83.80.Jx, 47.20.Ma 



I. INTRODUCTION 

The classic instability of hydrodynamic fingering oc- 
curs when a fluid of a certain viscosity is injected into a 
more viscous fluid between closely-spaced parallel plates 
[l|; this Saffman- Taylor instability Q has been widely 
studied, and has had many variants, including fingering 
in polymer fluids, foams and gels [1, 0, H, E 0|- The 
morphology of the instability can also be influenced by 
anisotropy Q , wetting Q or modifications of the surface 
tension or densit y co ntrast between fluids due to a chem- 
ical reaction [IfllTlT] . Hydrodynamic Angering is one of a 
wider class of instabilities occurring when one material is 
injected into another, or grows from a chemical or biolog- 
ical process; other such systems include smoldering flame 
fronts [lH , filamentary microorganisms [l3| , silica or iron 
tubes forming around metal salt solutions [3, , elec- 
trochemical deposition of metals growing biofilms 
[IS EH) or l ava nows 0- I n many of these cases, the 
interface itself is defined by a reaction or a solidification 
which changes the material properties, and is driven by 
the expanding growth of the interior. The patterns in 
these systems result from a complex interaction of reac- 
tion, diffusion and mechanical effects in association with 
a particular rheology or elasticity of the interfacial zone. 

Here we present a new type of pattern-forming insta- 
bility, where fingering is linked to a rheological change - 
a gelling - due to a reaction which produces a viscoelas- 
tic micellar medium at the interface between two aque- 
ous solutions of identical viscosities in a Hclc-Shaw cell 
(20| . The micellar gel-like medium is formed by reaction- 
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diffusion at the interface between two otherwise ordinary, 
miscible water-like Newtonian fluids of identical viscosi- 
ties. It is worth mentioning that we use the word inter- 
face in a more general sense than as immiscible, sharp 
boundary. Here, the interface is a region which thick- 
ens in time and shows no discontinuity on a macroscopic 
scale. One distinguishing aspect of our system is that 
the fingering occurs independent of which fluid is injected 
into which, making this a true interfacial instability. In 
this paper we present our first exploration of this insta- 
bility by varying fluid characteristics (concentrations) as 
well as the flow rate. 



The fluids we study are aqueous solutions of the 
cationic surfactant cetyltrimethylammonium bromide 
(CTAB), and the organic salt sodium salycilate (NaSal), 
which form a strongly viscoelastic micellar material when 
brought into contact. It is well known that a viscoelas- 
tic fluid is produced by the assembly of surfactants into 
long wormlike micelles, driven at low volume fractions 
by the mediating presence of certain slightly hydropho- 
bic organic counterions [U [22j • The essential difference 
between these wormlike micellar fluids and more stan- 
dard polymer fluids is that these aggregates are in a dy- 
namic equilibrium with free surfactants in solution due 
to breaking and reforming processes (23j . 



Much of the work on wormlike micellar fluids has fo- 
cused on their unusual properties, either Theologically 
(USE HE S3 or hydrodynamically [H, HJ H [M[ ; there 
has been less study of the reaction process between the 
surfactant and the salt which gives rise to the wormlike 
micelles [Hj]. In fact, the simple act of preparing these 
fluids presents a startling phenomenon: two dilute so- 
lutions with essentially the physical properties of water 
become, upon combination, a strongly elastic fluid. 
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II. EXPERIMENTAL SETUP 

The experiment consists of injecting one of the solu- 
tions into a Hele-Shaw cell previously filled with the other 
solution. Two cells where used in two different laborato- 
ries. Cell 1 (PSU) is made of two 9.5 mm-thick square 
glass plates, 15 cm x 15 cm. Cell 2 (UJF) is made of 
two 10 mm thick glass disks, 20 cm in diameter. A gap 
b = 0.6 mm between the two plates is fixed by brass 
shims. A light box below the cell provided a nearly uni- 
form illumination (light box from Schott-Fostec for cell 
1, a4 electro-luminescent panel from Selectronic for cell 
2), and the pattern growth is observed from above by 
a COHU 4912 video camera connected to a Macintosh 
computer. 

In most experiments, the cell is initially filled with a 
solution of NaSal ( Sigma- Aldrich) of a given concentra- 
tion, after thorough cleaning of the glass plates so that 
the system is completely prewetted by the solution with 
no bubbles or wetting defects. A small amount of blue 
(McCormick) or green (Vahine) food coloring is added to 
this solution for purposes of visualization. The injected 
surfactant solution (CTAB, Sigma- Aldrich) is uncolored. 
In all experiments, both solutions have the same concen- 
tration, ranging from 25 to 50 mM. Injection is made 
through plastic tubing at the center of the top plate of 
the Hele-Shaw cell by a KDS100 syringe pump at a con- 
stant flow rate Q, which ranges from to 1000 ml/h 
(0.278 ml/s). In some experiments, the opposite arrange- 
ment of fluids was tested: NaSal solution injected into 
CTAB. 

Two kinds of experiments were run: i) radial injection 
experiments, where the fluid is injected in an axisym- 
metric fashion in the Hele-Shaw cell, and ii) linear ex- 
periments where the initial direction of flow was imposed 
thanks to a U-shaped shim placed around the injection 
hole, aimed at studying isolated fingers. 



III. EXPERIMENTAL RESULTS 

A. Radial geometry 

The two solutions of CTAB and NaSal, at equimolar 
concentrations ranging from 25 to 50 mM are completely 
Newtonian before they are brought into contact with each 
other, each with a viscosity indistinguishable from that of 
pure water. Consequently, in the absence of the "worm- 
like micellar reaction" at the interface, a Saffman- Taylor 
experiment would exhibit a completely stable radially ex- 
panding circle. Instead, as the front moves away from the 
injection point and the injected volume increases, a se- 
quence of growth regimes takes place. In addition, both 
fluids are miscible with no surface tension, and no wet- 
ting force exists at the moving front. One can therefore 
neglect wettability effects that play an important role in 
the Saffman- Taylor experiment with immiscible fluids Q . 




FIG. 1: The onset of the fingering instability for an expanding 
circular front driven by an injection rate of 150 ml/hr, 30 mM 
CTAB solution injected into 30 mM NaSal solution: t = 10.0 
s, t = 20.0 s, t = 28.8 s, t = 64.7 s. Scale: picture width is 
8 cm. 




FIG. 2: Fingering instability under the same conditions as 
Fig. [T] but with the opposite arrangement of fluids - 30 mM 
NaSal solution injected into 30 mM CTAB solution: t = 8.2 
s, t = 11.0 s, t = 16.3 s, t = 54.8 s. Scale: picture width is 
8 cm 



We begin by discussing the injection of the surfac- 
tant solution (CTAB) into the solution with the organic 
counter- ion (NaSal). A stable circular expansion is ob- 
served in our experiment at high injection speeds, which 
occur close to the injection point due to the constant 
flow-rate condition. In that case, when the radius R of 
the expanding CTAB solution is small, or equivalently 
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FIG. 3: Fingering patterns for different flow rates after the 
same total volume (4 ml) of 50mM CTAB has been injected 
into 50 mM NaSal: a) 20 ml/h; b) 50 ml/h; c) 100 ml/h; d) 
200 ml/h. The dark borders are the gelled interfaces, which 
bound the interior, flowing conduit of fresh fluid which drives 
the further growth of the pattern. Scale: picture width is 
11 cm. 

when the speed dR/dt is large, the gelling reaction oc- 
curs at a slower timescale than the stretch rate of the 
gelled membrane which remains thin and offers no resis- 
tance to flow: the boundary between fluids is stable and 
isotropic. The gel at the interface is probably very thin 
and dilute, and does not break or split (Fig.QJi). 

When R reaches some critical value which depends 
on flow-rate and concentrations, perturbations are ob- 
served on the smooth front (Fig. [TJa). The result of 
this instability is that the constant flow rate Q focuses 
into the perturbations, which become more "active" and 
bulge out into mushroom shapes. Meanwhile, the regions 
between the bulges slow down, and appear to harden. 
Subsequently, each individual mushroom spreads almost 
isotropically, albeit at a slower rate than the initial 
front, as each is effectively fed with a reduced flow rate 
(e.g. about Q/10 in Fig. QJ;). The appearance of these 
mushroom shapes around the perimeter bear a striking 
resemblance to the "breakout" of a lava flow [l9[ . 

A precise determination of the onset of the first in- 
stability leading to the mushroom pattern can be made 
by measuring the radius of the circle circumscribing 
the growing front as a function of time; when the cir- 
cular front becomes unstable to fingerlike protrusions, 
the growth rate of the gel front (and the entire pat- 
tern) goes through a rapid increase. We find an initial 
quantitative agreement with the constant flow rate law 
R(t) = y/ Qt/irb, where Q is the imposed flow rate and 
b is the gap width. As the gel front slows, an instability 
appears nearly simultaneously around the circular front. 
These perturbations grow, and take on the appearance of 
mushrooms, as shown in Fig. [T] These new fronts, which 
have focused the flow from the center, move at a faster 
rate than the front just before the instability. Eventually 
these fronts slow down as the pattern expands in size, and 
a second generation of fingers breaks from the pattern. 

Surprisingly, we observe the same instability to the 
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FIG. 4: Phase diagram for the injection of 50 mM CTAB 
into 50 mM NaSal: Time at which successive regimes are 
triggered vs. flow rate. (•): onset of mushroom growth, (A): 
onset of fans, (o): onset of tentacle growth. Dashed lines are 
guidelines. 



mushroom pattern by performing the inverse experiment: 
injecting a 30 mM NaSal solution at the same rate into 
a 30 mM CTAB solution (Figure [2}. The fact that the 
instability occurs independent of which fluid is injected 
into which distinguishes this from Saffman- Taylor and 
other related fingering instabilities, and shows that this 
fingering is truly driven by an interfacial instability. Nev- 
ertheless, the critical radius at which the circular front 
becomes unstable is significantly smaller as can be seen 
when comparing Figures [1] and O This may be due to 
different diffusion coefficients of CTAB and Salicylate, or 
different compositions of the micellar media formed when 
injecting the surfactant into the salt or the opposite ar- 
rangement. 

In the mushroom instability, the interfacial membrane 
between the two solutions appears to break simultane- 
ously at a number of points along the perimeter. Bulges 
of fluid protrude, and the gel-like material forms at the 
new interfaces. When this becomes more solid, it resists 
further smooth growth, but under the imposed condi- 
tions of constant injection rate, the pressure of the inter- 
nal fluid will increase until a rupture of the gel membrane 
occurs again, and the process of fingering repeats. 

Close observation of the interface between the two flu- 
ids reveals that the membrane is thickening with time, 
evidently due to the further reaction of surfactant and 
organic counterion diffusing through the material already 
formed. This thickening is made obvious by the slightly 
darker color of the material at the interface. It has also 
been confirmed by a 'post-mortem' test after the exper- 
iment: when the glass plates are separated, that the 
darker material is indeed gel-like, and clings to the glass 
while the two fluids rapidly flow away as the cell is emp- 



FIG. 5: A regular pattern of fingers in the fan state, obtained 
when injecting 35 mM CTAB into 35 mM NaSal, Q = 600 
ml/hr. Scale: picture width is 10 cm. 

tied. This is most likely the reason that the membrane 
elasticity and resistance to flow increase with time. When 
some material threshold is reached [33[ , a fracture occurs 
in the membrane and fresh fluid leaks and forms the next 
generation of mushrooms. However, the relatively well- 
defined wavelength of the mushroom pattern and the si- 
multaneity of finger appearance suggests an underlying 
mechanism based on gel characteristics (thickness, con- 
centration) prior to fracture. 

Note that the subsequent evolution of the patterns in 
Figs. 1 and 2 differs for the two injection arrangements: 
while the experiment injecting CTAB solution into NaSal 
solution leads to other instabilities (described next), the 
other arrangement is simpler. The same process of ex- 
pansion, hardening, fracture, and break-out seems to re- 
peat as the pattern expands, without any other mor- 
phologies. 

We have identified two other instabilities which occur 
at later times as the radial pattern expands. After a few 
generations of mushroom patterns, we observe a tran- 
sition to a fan-like pattern of contiguous wide fingers, 
giving the growing front a sort of flower-like appearance. 
The fingers emerging from these fans eventually become 
narrower "tentacles", which then undergo a curling in- 
stability and start to meander. These tentacles seem to 
prefer growing along an existing structure rather than 
into fresh fluid. By changing the injection rate Q, the 
different growth regimes can be shifted both radially and 
in time: at higher injection rate, a given pattern instabil- 
ity is found to occur at larger distances from the center 
while at lower flow rates, the isotropic regime is almost in- 
visible and the first observed pattern can be mushrooms 
or even multiple tentacles diverging from the injection 
point (see fig. [3]). The way this sequence of instabili- 
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FIG. 6: Critical radius of breakdown of axial symmetry vs. 
flow rate for three different CTAB / NaSal concentrations. 



ties and patterns takes place as a function of flow rate is 
summarized in the schematic phase diagram of fig. 2] for 
50mM concentrations. 

While some regimes can be skipped by lowering the 
flow rate, which effectively shrinks patterns radially down 
to the center, there are also ranges of concentrations and 
flow rates for which the fan instability is the first one 
to occur, directly after a stable circular front has grown 
to a finite distance to the center, as shown in Fig. [5] for 
Q = 600 ml/h and 35 mM concentration. 

The second pattern observed in this experiment 
( "fans" ) clearly involves an instability with a well de- 
fined wavelength (figs. [3Jd) andO; while it appears to 
be similar to the first instability, the wavelength is much 
smaller. This underlines a major influence of the radial 
expansion: a consequence of interface elongation is an az- 
imuthal tension in the gel membrane which stabilizes the 
shortest wavelengths, in a similar fashion to surface ten- 
sion in the classic Saffman- Taylor instability As the 
radius of the injected fluid increases, membrane stretch 
rate decreases and the interface destabilizes. 

More quantitative data on the onset of the first insta- 
bility seen after the stable circular growth was obtained 
from the conditions when the stable circular growth stops 
and either a regular array of fingers characteristic of the 
fan state (as in Fig. |5}, or a mushroom pattern (as in 
Fig. [T|) appears. Figure [5] shows the evolution of the 
critical radius at which the axial symmetry is lost vs. 
flow rate for three different concentrations. The critical 
radius increases with the flow rate: the gel membrane 
stretch rate being higher, it takes longer to become thick 
enough and trigger the instability. As far as concentra- 
tion is concerned, one would intuitively expect that the 
stronger the concentration, the harder the gel. Thus the 
gel breaking should occur sooner for higher concentra- 
tions. However, an inversion is seen in Fig. [SJ the curve 
for 25 mM is below others at high flow rate. This might 
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FIG. 7: Growth of the micellar gel pattern in the tentacle 
state, for the injection of a 50 mM CTAB solution into a 50 
mM NaSal solution, at a flow rate Q = 20 ml/h. The width 
of each image is 11 cm. 



be due to a different rheology at different concentrations 
(e.g. shear thinning occurring at different shear rates). 

The evolution of fingers in the tentacle state (Fig. [7]) 
is significantly different than either the mushroom or fan 
state; the growth is very localized at the tip. The fresh- 
est fluid arrive at the tip from the interior of the finger, 
where it has not yet hardened to an immobile state, and 
is advected to the sides by the feeding flow, where the 
interface solidifies. We expect this instability to occur 
when the rate of creation of new interface exceeds the 
hardening time of the gel. 

A remarkable feature of these tentacles is their roughly 
constant width. Since in this regime the average number 
of actively growing tentacles Nf is typically constant (see 
Fig. [SJ which appears to have reached a steady state of 
~ 11.4 fingers), the velocity of the moving tip must on 
average also be constant. Tentacles can therefore be seen 
as a state to which the system is attracted under certain 
conditions are reached: critical width (or wavelength of 
the "fan" instability) and critical velocity. 

The stability of this regime can be qualitatively ex- 
plained in the following way: a slight decrease of the fin- 
ger width leads to a velocity increase because of the con- 
stant flow rate. The gel at the tip would then be weaker, 
which favors radial growth and tip widening. Conversely, 
a slight increase in finger width would decrease the fin- 
ger velocity, which would give more time for the gel to 
harden. This would result in a greater resistance to fin- 
ger widening. Gel fracture would then occur and redirect 
the flow. We focus on this finger regime in next section. 

Curling structures appearing in the ultimate develop- 
ments of tentacle growth, the tip favors motion along an 
existing membrane. One mechanism for this could be 
that in these regions a depletion of the outer solution 
(NaSal) has occurred. As a consequence, the gel forming 
there is weaker and poses less resistance to flow, resulting 
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FIG. 8: The total number of actively growing fingers Nf in 
the tentacle growth state as a function of time (50 mM CTAB 
solution injected into 50 mM NaSal solution at 20 ml/hr). 



in an effective attraction between growing fingers. 



B. Linear geometry - single finger growth 

While our initial observation of micellar interface fin- 
gering was made in a radial injection experiment, the 
geometry complicates the dynamics at a fixed flow rate, 
since the front speed decreases with radius as the pattern 
spreads out. We therefore performed a second series of 
experiments in a linear geometry at LSP in Grenoble. 

In order to investigate the remarkable finger growth 
regimes seen in the ultimate stages of radial growth, we 
made experiments focused on the study of this particular 
pattern. A small device (U-shaped plate) was placed in 
the Hele-Shaw cell around the injection point to promote 
single finger growth in a certain direction, allowing pre- 
cise control of the flow rate inside the finger. The plate 
forms a small channel of width 3 mm and length about 
2 cm. Note that this channel serves only to set the ini- 
tial condition: finger growth is observed once the fluid is 
coming out of the channel: the finger grows in an uncon- 
fined half-space. When the flow rate Q is in the correct 
range (approximately 0-3 ml/h), a single finger grows. 
For Q > 3 ml/hr, it becomes unstable and splits into two 
or more fingers. 

We observed different types of fingers, depending on 
flow rate and concentration see Fig. [9l For a given con- 
centration, at low flow rates fingers are evanescent: their 
width decreases as the finger grows until the narrow tip 
opposes too much resistance to flow. Then the membrane 
breaks somewhere and a new evanescent finger grows. At 
intermediate flow rates (close to 2 ml/h), there exists a 
narrow range where steady growth occurs: fingers keep 
a constant width and the front moves at constant speed. 
Above this range, one observes surprising oscillating fin- 
gers: the front periodically expands and contracts, mod- 
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FIG. 9: Linear finger growth regimes. Left to right: evanes- 
cent finger (concentration c = 30 mM, flow rate Q — 1 ml/h), 
steady finger (c = 40 mM, Q = 2.5 ml/h), wavy /oscillating 
finger growth (c = 30 mM, Q = 1.5 ml/h). Each picture is 2 
cm wide, and the time interval between wavy finger pictures 
is 150 s. 



ulating tip width. 

Note that in both cases only the finger tip is active; the 
sides remain fixed, and only thicken slightly with time as 
discussed above. Another remarkable feature is the tip 
shape: it is a straight front, perpendicular to the direc- 
tion of motion, except for evanescent fingers at the end 
of their life. This is completely unlike the rounded tips 
of most Saffman- Taylor fingers, and also the needle-like 
fingers seen in fingering experiments on non-Newtonian 
associating polymer fluids [34] or colloidal pastes [3o| . 

Moreover, the fact that the finger has a selected width 
while being far from any sidewalls is very different than 
the classic Saffman- Taylor case or the fan state described 
above, where the finger width is half of the channel width 
(without surface tension effects) [l| ; in a radial flow with 
far away boundaries the fingers expand and split [361 [37| . 
This difference is a consequence of the rapid formation 
of a gel-like membrane that poses a strong resistance 
to flow, preventing finger widening. We suspect that in 
evanescent fingers the pressure and velocity are too low 
to prevent quick widening of the gel membrane toward 
the inside of the finger. The oscillation mechanism in 
wavy fingers is more mysterious. The oscillation mecha- 
nism in wavy fingers is more mysterious. Given the fact 
that timescales involved in this process (~ 100 s) could 
be of the same order of magnitude as rheological times of 
the micellar fluid, it may share a common origin with the 
stick-slip instability suggested by Puff et al. [a] , however 
a more detailed characterization of the material will be 
needed to test the connection. 



IV. DISCUSSION AND CONCLUSIONS 

We have presented here a new kind of fingering insta- 
bility in a Hclc-Shaw cell, defined by the interface be- 
tween two solutions which together form wormlike mi- 



celles. The most striking aspect of these fingers is that 
the instability is completely determined by the properties 
of the interface. This is shown by the fact that essentially 
the same instability occurs independent of which solu- 
tion is injected into which (in fact both solutions have 
the same viscosity). We have seen that there are some 
differences between these two patterns, which may arise 
from an asymmetry in diffusion coefficients of the two re- 
acting species - it is very likely that the salicylate ion has 
a higher mobility in water than either a spherical micelle 
or a single surfactant molecule. 

To our knowledge, the complex morphology of a grow- 
ing, gelling interface has not previously been studied; 
however the early work by Hatschek on the gelling of 
a sinking drop of reactive fluid during vortex formation 
is similar in spirit [38l . |39| . From our observations, many 
physical processes seem to influence pattern formation 
in this system: reaction-diffusion, flow and rheology of 
the viscoelastic micellar gel. The instabilities of such a 
system provide a new set of mechanisms which are of po- 
tential relevance to processes occurring in the growth of 
biological structures like bacterial biofilms. 

Although our system is not explicitly biological, it in- 
cludes many physical aspects of biofilms, such as the for- 
mation and growth of an elastic gel subjected to flow. 
Our study sheds light on the purely physical instabilities 
to which such a system is susceptible. The general prob- 
lem of an expanding or growing gel finds an important 
application in certain biomedical infections, where the 
bacteria secrete a resistant "biofilm" matrix, which sup- 
ports and protects them, often against disinfectants and 
antibiotics; an important example is Pseudmonas Aerug- 
inosa, responsible for many hospital infections [40]. Our 
experimental results seem most closely related to other 
finger-like structures driven by dynamic growth or cool- 
ing processes, such as those produced by the rupture of 
the cooled crust of flowing lava [l9(. While much re- 
mains to be done to reach a quantitative understanding, 
the experiment we present here shares common qualita- 
tive features with many of these systems, and could be 
seen as a representative reaction-diffusion-advection sys- 
tem, displaying similar pattern formation but in a more 
convenient to controlled laboratory study than infectious 
biofilms or flowing lava. 
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